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Abstract 



O 

O , 

C^l , Quantum computers are considered as a part of the family of the re- 

versible, lineary-extended, dynamical systems (Quanputers). For clas- 
■ sical problems an operational reformulation is given. A universal algo- 

rithm for the solving of classical and quantum problems on quanputers 
£S) \ is formulated. 

^ ■ 1. Time-reversible classical discrete dynamical systems (comput- 

ers) and corresponding (quantum) extraparts 

The contemporary digital computer and its logical elements can be consid- 
ered as a spatial type of discrete dynamical systems: 

(N 

o; s n (k + i) = (j> n (s(k)), (l) 

Qh! where s n (k), 1 < n < N(k), is the state vector of the system at the discrete 

time step k. Note that, with time steps k, can change not only value, but also 
ctf the dimension N(k), of the state vector s(k). 

^ { - Definition. We assume that the system ([!]) is time-reversible if we can 

define the reverse dynamical system 



s n {k) = <p-\s{k + l)). (2) 
In this case the following matrix 

_ dMs(k)) 

nm ~ ds m (k) ' (3) 

is regular, i.e. has an inverse. If the matrix is not regular, this is the case, for 
example, when N(k + 1) ^ N(k), we have an irreversible dynamical system 
(usual digital computers and/or corresponding irreversible gates). 

Let us consider an extension of the dynamical system ([!]) given by the 
following action functional: 

A = J2Uk)(s n (k + l)-Ms(k))) (4) 

kn 



1 The results of this paper where presented on the Workshop 
"Quantum physics and communication ", Dubna, May 16-17, 2002. 
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and corresponding motion equations 



dH 

s n (k + 1) = (f) n {s{k)) 



dl n (k) ' 

l n (k - 1) = l m (k) ^0^ = l m (k)M mn (s(k)) = J^-, (5) 

where H = J2kn ln(k)(j) n (s(k)), is discrete Hamiltonian. In the regular case, we 
put this system in an explicit form 

s n (k + 1) = <f> n (s(k)), 

l n (k + l)=l m (k)M-l l (s(k + l)). (6) 

From this system it is obvious that, when the initial value l n {ko) is given, 
the evolution of vector l(k) is defined by evolution of the state vector s(k). So 
we have the following 

Theorem. The regular dynamical system ([I]) and the extended system (0) 
are equivalent. 

Note that the corresponding statement from M concerning to the continual 
time dynamical systems is not so transparent as our statement for discrete time 
dynamical systems. In the continual time approximation, the discrete system 
(H) reduces to a corresponding continual one, @: 

dv 

Xn{t) = V n (x), p n (t) = —Q^rPm- (7) 

Indeed, let us change the dependent variables, 

s n (k) = x n (t k ), l n (k) = p n {t k ), t k = kAt, At « 1. (8) 
Then the action functional @ can be transformed as follows: 

A = Y,Pn(t k )(x n {t k + At) - <j) n (x(t k ))) 

kn 

Atp n (t k )(x n (t k ) - V n (x(t k ))) 

kn 

dtp n (t)(x n - v n (x)), v n (t k ) = (0„(x(4)) - x n (t k ))/At, (9) 



and corresponding motion equations are presented by the system (|7|). The 
equation of motion for l n (k) is linear and has an important property that a 
linear superpositions of the solutions are also solutions. 

Statement. Any time-reversible dynamical system (e.g. a time-reversible 
computer) can be extended by a corresponding linear dynamical system (quantum- 
like processor) which is controlled by the dynamical system and which, due 
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to the superposition and entanglement properties, has a huge computational 
power. 

Note that || for the following (infinite - dimensional) equation of motion 

iV t = AV- l -V\ (10) 

the corresponding linear subsystem is nothing but the familiar Schrodinger 
equation, 

iipt = -Aif) + Vip. (11) 

The extended system 

iV t = AV- ^V 2 } 

iip t = -Ail) + Vip, (12) 

is mathematically and physically complete, as potential is not an external 
function, but is part of the Hamiltonian dynamical model. 

2. Solution of the Quantum problems with the Methods of the 
Quantum Computing 

The standard classical computing - the technology of storing and transform- 
ing information, is based on the classical physical theory and can be universally 
divided (factored) on the hardware (body) - memory and processor of comput- 
ers, and software (soul) - algorithms and programs of users, parts. Under the 
Quantum computing (Quanputing, QC), we mean an extension of the classical 
computing by hardware which is described by the (corresponding extention of 
the) quantum theory - quantum registers and quantum logical units and corre- 
sponding software. In Quanputing the factorization on hardware and software 
is only an approximate concept, < SH >=< S >< H > +0(h). 

A quantum system can be described by corresponding Schrodinger equation 

I 

ih^\i/>>=H\1>>, (13) 

where \ip >= \ip(t) > is the state vector from the state Hilbert space and 
H = H(p, x) is an operator-Hamiltonian. In the case of one nonrelativistic 
particle the operator isH = T + V = j^ + V(x). The fundamental bracket is 
[x,p] = xp — px = ih. The configuration space form of Eq. ( |13"D is 

ih^^ = H^(x,t), (14) 

where 4>(x,t) =< x\ip(t) > and Hamiltonian is H = — + V(x). 
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The formal solution of Eq. (|13|) is 

|V(f)>=l7(*)|Vo>, U(t)=exp{- % -tH). (15) 

n 

The main steps maid in QC are the following: 

U{t) = {U X / N ) N = (U T U V ) N + 0{1/N), (16) 

where 

U l ' N = exp {-OH) = exp (-6T)exp (-6V) + 0(1/N 2 ) 

= U T U V + 0(1/N 2 ), e= % -r, T = ± (17) 

Then, for a corresponding matrix element we have (see Appendix) 

< x n+l \U l ' N \x n >~ exp (g(^( Xn+1 ~ Xn ) 2 - V{x n ))) + 0(1/N 2 ) (18) 

and 

% ^ Ttl X X 

< x out \U(t)\x in > ~ / dx 1 dx 2 ...dx N exp (-r V( — (-^ -) 2 - V(a; n ))) 

+ 0(1/N), \x >= \x in >, < arjv+il =< Xoatl- (19) 

This finite dimensional integral representation of the matrix element is in 
the ground of the functional (continual) integral formulation of the quantum 
theory 0]. Then we discretize the wave function in flnp, ([TJJ) and impose the 
periodic boundary conditions || a n (t) = if)(x n ,t), a n+ ^ = a n , x n = nAx. 

We store these amplitudes in a k-bit quantum register, 

N 

|V(t) >= £a n |n>, N = 2 k , (20) 

n=l 

where \n > is the basis state corresponding to the binary representation of 
the number n. The second factor in (fP7D in a coordinate representation corre- 
sponds to a diagonal unitary transformation of the quantum computer state 
if)(x,t). After Fourier transformation ip(x,t) into momentum-space represen- 
tation ip{p, t) the first factor in flT7| ) can be applied in the same way. Diagonal 
unitary transformations of the type |n >— > e jF(ra )|n >, where F(n) is a function 
of n, can be done || with the following succession of steps 

\n> -> |n,0 >^ |n,F(n) >-»• e iF(n) K^(™) e iF{n) |n,0 > 

-> e iF(n) |n>- (21) 
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For n particles in d dimensions we need nd quantum registers. If the po- 
tential V(xi, 22, ■••) couples different degrees of freedom, we need the diagonal 
unitary transformations acting on several registers, 

\n u n 2 , ... >^ e^'^'-Vi, n 2 , ... > . (22) 

The disctetized bosonic quantum fields are naturally defined on the space- 
time lattice, scalar (particles) fields - on lattice sites, vector particles - on 
lattice links, tensor particles - on lattice plackets, and so on (see e.g. 0). 

3. Solution of the Classical Problems with the Methods of the 
Quantum Computing 

If we can solve the quantum problems with the quantum computers ||, 
we can probably solve the classical problems, too. With the paper || we 
started the investigation of computational hard problems of classical physics, 
such as the turbulent phase of hydrodynamics, with the methods of Quantum 
Computing (QC). 

We usually formulate the classical problems as a Hamiltonian system of 
motion equations [|TIJ 

dH dH 

<in = TT- , Pn = - T -, 1<U<N, (23) 

where q n are coordinates for the configuration space of the system and p n are 
corresponding momentum, H = H(q,p) is a Hamiltonian function(al). The 
system (^3]) belongs to a more general class of the dynamical systems defined 
by the following system of motion equations 

x n = v n (x), 1 < n < N, (24) 

when N is an even integer and 

_ dH(x) 

OXvry 



'Tim- 



with the fundamental canonical bracket {x n ,x m } = e r 

Note that any dynamical system ( p4] ) can be extended to the following 
Hamiltonian form, (see, e.g. and section 1 of this paper.) 

X n V n [x) (in) Hfi 

Ov 

Pn = '-Q^rVrn = {Pn, H} , (25) 

with Hamiltonian H (x, p) = J2 n v n{x)Pn and the fundamental canonical bracket 

For any observable i/j(x) of the Hamiltonian dynamical system, we have the 
following motion equation: 

4j(x) = ty, Ht} = enm^^ 1 = ~U = ~H 2 i>, (26) 
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where 



Vn ^~dx~ n = n Vn ^ n = li H2) ^ 



Hi is the classical Hamiltonian function and H 2 is the quantum Hamiltonian 
operator. Now we can consider Eg . (p6|) as a classical Hamiltonian one with 
Hamiltonian Hi or as a quantum Schrddinger equation with the Hamiltonian 
operator H 2 . For the second form we can apply the QC methods developed in 
the previous section. 

Note that the operator H 2 is a self-adjoint operator, H 2 = H 2 , when 
divv = and we can put it in the second-order form with respect to the 
momentum operator, 

1 i 
H 2 = v k p k = -(pkV k + v k pk) = -[pi, u k ], (28) 



2 V " r ' fx 



where 



Xu 

1 



Uk = 2 J ^ XkVk ' ^ — ^ — ^' ^ — 2- (29) 
The formal solution (|15|) is 

i - t 
\i){t) >= U(t)\ifo >= exp (--tH 2 )\^ >= exp (-- [u,p 2 })\ip > 

n n 

and for the factorized form ([16]) we have 

d 

U^ N = exp (-r 2 [u,p 2 ]) = J] e'^e^e-^e"^ + 0(r 3 ), (30) 

k=i 

where 

y h 2 N' 

Now we are ready to apply the formalism of the QC considered in Sec. 3 
to the classical problems. In the extended version of this paper, we describe 
some computational problems from different parts of physics. The Hamiltonian 
dynamics of heavy-particle accelerators; the standard model of condensed state 
physics (Hubbard model); the 123 standard model of high-energy (elementary 
particles) physics; unified string and M-theory models. 

4. Conclusions and perspectives 

In this paper we constructed an algorithm for solving difficult computa- 
tional(hard) problems of classical and quantum physics on quanputers. Note 
that before we get a real quanputer, we can simulate it on classical computers 



1% . It is interesting to investigate the functional integral formulation of the 
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classical theory based on a discrete representation (|30|). Then, Dirac's equation 
for electron 

ihd t i)(t,x) = (a n p n + (3m)ip(t,x) = {=\pl,u n ] + fimjip, (31) 

(where u n — | / dx n a n = a n x n /2,) for m — 0, has the form similar to the 

classical equation (|26D , (with a n as velocity matrices) and the Planck's constant 
h can be cancelled. This is an inverse to the consideration made when we try to 
find quantum (photon) nature of the (classical) Maxwell's equations and put 
them in the "Dirac's form", (see e.g. [0). It is possible to formulate similar 
equations for high spin massless fields. It is interesting to find a corresponding 
classical interpretation of these equations and make observable predictions. 
The effects of two of them (the photon's and graviton's) are obvious in everyday 
life. What about the effects of other massless particles? 

Appendix 

1. In the main text of this paper we used the following relation: 

e sA e eB e -sA e -eB = ^[A,B] + Q ^y ^ 

The relation 

(1 + eA)(l + eB)(l + eA)-\l + eB)- 1 = (1 + e 2 [A, B}) + 0(e 3 ), (33) 

may also be useful. 

2. Coordinate and momentum state vectors are correspondingly \x > and 

\p>, x\x >= x\x >, p\p >= p\p >, 

1 i d 

< p\x >= ip x {p) = ■— - exp {--px), xip x {p) = ih—ij x (p) = xij) x {p), 

vZnh dp 

< x\y >= [ dp < x\p >< p\y >= [ exp (^p(x — y)) 

J 1 J znh h 

= 5(x-y). (34) 
Now we calculate the following matrix element 

< %n+i\ exp (— ap 2 )\x n > = J d D p < x n+ i\p >< p\x n > exp {—ap 2 ) 

/dp . .p(x n +i x n ) o\ A , (x n _|_i x n ) . / r\ 

j2^r expil — n ap) = j2^r expi — 4^^ } ' (35) 

where in the case of the quantum mechanics of the particle, (|lj 



and in the case of classical mechanics, (|3 



a = t (J-) 1 /2 (37) 



A = J dp exp (—ap 2 ) = ^J—. (38) 
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